ELLIPTIC SYSTEMS AND MATERIAL INTERPENETRATION 



GIOVANNI ALESSANDRINI AND VINCENZO NESI 

Abstract. We classify the second order, linear, two by two systems for which 
the two fundamental theorems for planar harmonic mappings, the Rado- 
Kneser— Choquet Theorem and the H. Lewy Theorem, hold. They are those 
which, up to a linear change of variable, can be written in diagonal form with 
the same operator on both diagonal blocks. In particular, we prove that the 
aforementioned Theorems cannot be extended to solutions of either the Lame 
system of elasticity, or of elliptic systems in diagonal form, even with just 
slightly different operators for the two components. 



1. Introduction 

"A basic requirement of continuum mechanics is that intcrpenetration of matter 
does not occur, i.e. that in any deformed configuration the mapping giving the 
position u{x) of a particle in terms of its position x in the reference configuration 
be invertible". J. M. Ball 0. 

It is then a natural question to ask what systems of equations, among those 
used as models for elastostatics, give rise to invertible solutions when reasonable 
boundary conditions are prescribed. 

In this note we shall prove by an example that the Lame system of isotropic, 
linearized elasticity in the plane, with constant Lame coefficients, may lead to 
physically unacceptable solutions, because interpenetration of matter occurs. Let 
us recall here that the same phenomenon was previously found by Fosdick and 
Royer-Carfagni |12j for a more involved anisotropic linear system, by elaborating 
on an example due to Lekhnitskii [13] - In higher dimensions similar phenomena 
occur. From the mathematical point of view a basic example is due to De Giorgi 
[10 . In all these examples, however, a basic common feature is the presence of 
some sort of point singularity in the solution itself (in dimension greater than two) 
or at least in the gradient in any dimension. Such a singularity can only be present 
when the coefficients of the elliptic system are irregular, for instance discontinuous. 
Our examples are different from those previously known in several ways, but the 
most crucial difference is that we choose smooth and in fact constant coefficients 
in our systems. Therefore our examples, besides bringing a new argument to the 
many already well known, see for instance Ciarlet [8, p. 286] and the references 
therein, about the limitations of linearized elasticity, shed some new light on the 
tightness of certain classical properties enjoyed by harmonic mapping showing that 
they cannot be easily extended even within the class of constant coefficient systems. 

We now recall some fundamental properties of planar harmonic mappings, see 
the book of Duren [TT] for a very broad treatment of this subject. 
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We begin with the classical theorem of Rado, Kneser and Choquet. This theorem 
which was first stated by Rado [22], proved by Kneser [TB] immediately after and 
then independently rediscovered by Choquet [7J , remains the basic unequalled result 
of invertibility for mappings solving an elliptic system of equations. 

Let B be an open disk in the plane, let $ : dB — > 7 C M 2 be a homeomorphism 
of dB onto a simple closed curve 7. Let u € C 2 (B,R 2 ) U C(B,R 2 ) be the solution 
to the Dirichlet problem 

, . ( Au = , in B, 

1 1 \ m = $ , on dB. 

Let D be the bounded region such that dD — 7. The Rado-Kneser-Choquet 
Theorem states the following. 

Theorem 1.1. If D is convex, then u is a homeomorphism of B onto D. 

The main reasons why this theorem remains substantially unequalled are 

(i) no analogue of this theorem holds true in dimension three or higher, as it was 
shown by a striking example by Laugesen [17], see also Melas [2T] . 

(ii) the convexity assumption on the target domain D is optimal. In fact, it is 
known since Choquet [7] , that if D is not convex then there exists homeomorphisms 
<f> : dB — > dD, for which the injectivity of u, the solution to fails. See also 
[2 J for a thorough investigation of this issue. 

Nevertheless various kinds of generalizations of the Rado-Kneser-Choquet The- 
orem have been obtained. Regarding harmonic mappings between manifolds, see 
Schoen and Yau |23j and Jost |15j . For mappings u whose components solve a linear 
elliptic equation, let us mention Bauman, Marini and Nesi [B] and also [T], [3J. It is 
worth pointing out that, in these last two papers the Rado-Kneser-Choquet The- 
orem has been extended to linear elliptic equations in divergence form with merely 
bounded measurable coefficients. For quasilinear equations of the p-Laplacian type 
see [4]. 

As a remarkable special case of our Theorem |1.4[ which will be stated here below, 
we prove that no analogue of the Rado-Kneser-Choquet Theorem holds when the 
diagonal Laplacian system is replaced by a Lame system with constant moduli of 
the following form 

H div((Vw) T + Vu) + A V(div u) = 0. 
More precisely, we have the following. 

Theorem 1.2. Let [i, A € K with \i > and n + A > 0. There exist a disk B C K 2 , 
a bounded convex domain Del 2 and a smooth diffeomorphism $ : dB — -> dD, so 
that the unique solution u £ W 1,2 (B, M 2 ) to 

f fi div((Vu) T + Vu) + A V(div u) = , in B, 
[ u — <f> , on dB, 

is not a homeomorphism of B onto D. 

Note that the natural unknown for the Lame system, should be the displacement 
field 5, rather than the deformation field u. However, due to the fact the the identity 
mapping / is also a solution of the Lame system, we have trivially that S solves the 
Lame system if and only u — I + 5 solves the same system. This is the reason why 
in Theorem |1.2| it is understood that the solution u is representing a deformation 
field. 
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In Remark |2.3| we shall see that the solution fails to be an homeomorphism in a 
very strong way, in fact u maps B onto a domain larger than D, moreover it folds 
itself along a curve, across which the orientation is reversed. 

Our next Theorem is a far reaching generalization of the previous one. It shows 
that the Rado-Kneser-Choquet Theorem holds only if one deals with elliptic sys- 
tems of diagonal form with the same scalar elliptic operator on both components. 
We need some definitions. Consider a constant coefficients second order elliptic 
system of the form 



(1.2) 



div(AVw 1 + BVu 2 ) = 0, 
div(CVw 1 + DVu 2 ) = 0. 



where A, B, C and D are 2x2 real constant matrices, and the unknowns u 1 and 



u 2 are real valued functions, we say that the system (1.2 1 is elliptic if it satisfies 
the Legendre-Hadamard condition 

(1.3) r£A£-£ + Thri2(.B + C!)€-€ + rftDt.t>0, for every ^,?)eR 2 \{0}. 

This condition is weaker than the strong convexity condition, namely the positivity 
of the 4x4 matrix given in block form as 



M 



A B 
C D 



and, as it is well known, ellipticity ( |1.3| is the same as rank one convexity of the 
quadratic form associated to M, see Theorem 5.3]. We note that there is no loss 
of generality in assuming that the matrices A, B, C and D are symmetric. 



Definition 1.3. We shall say that the system (1.2|) is equivalent to the system 

div(A'Vw 1 + B'Vu 2 ) ■■ 
div(C"Vw 1 + D'Vu 2 ) 

if 1 here exists a ni)ii->iii»iilar 2x2 matrix ( a ^ ) such thai 



0. 

= 



(1.4) 





B \ 




A c 


D J 


= ( 



aid (3Id 
7ld <5Id 



A' B' 
C D' 



Theorem 1.4. Let the ellipticity condition {1.3) be satisfied. The following alter- 
native holds. 
Either 



i) the system (1.2) is equivalent to the system 
(1.5) 



div^Vw 1 ) = 0, 
div(AVw 2 ) = 0, 



ii) there exist a disk Bel 2 , a bounded convex domain D CM 2 and a smooth dif- 
feomorphism $ : dB —>■ dD, so that the unique solution u = (it 1 ,?/ 2 ) € W 1 ' 2 (B, M 2 ) 
to 

div(AVu 1 + BVu 2 ) = , m B, 
div(CVu 1 + DVu 2 ) = , on B, 
u = <i> , on dB. 

is not a homeomorphism of B onto D. 
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The above results show that one of the most basic properties enjoyed by planar 
harmonic mappings cannot be extended to other elliptic systems in the plane. It is 
then natural to ask similar questions for another fundamental property of injective 
harmonic mapping. A benchmark of the theory is a result of H. Lewy |18j proving 
that harmonic homeomorphisms are, in fact, diffeomorphisms. More precisely we 
have the following result. 

Theorem 1.5. (H. Lewy.) Lut u — (u 1 ,^ 2 ) : B — > M 2 be a harmonic mapping. If 
u is invertible, then 

(1.6) det-Duy^O for every [x,y)£B. 

Also in this case the validity is limited to two dimensions. J. C. Wood [24] found 
a third degree polynomial harmonic mapping which provides a counterexample in 
dimension three. On the positive side, Hans Lewy [TS] recognized that, in three 
dimensions, if u is the gradient of an harmonic function and it is a homeomor- 
phism, then it is a diffeomorphism. This result was extended to any dimension in 
the remarkable paper by Gleason and Wolff 14J. In a different direction, several 
generalizations of Lewy's Theorem have been achieved in dimension two when the 
components of u satisfy the same linear elliptic equation of the form div((rVw J ) = 0. 
For the case of sufficiently smooth a see [B]. When a is allowed to be discontinuous, 
weak forms of Lewy's Theorem have been obtained in pQ and pj], A version for 
p-Laplacian type equations can be found in pQ. 

In the next theorem, we show, by means of examples, that Theorem [L5] cannot 
be extended to an arbitrary elliptic system with constant coefficients, unless, again, 



the systems has the special form (1.51 



Theorem 1.6. Let the ellipticity condition (1.3) be satisfied. The following alter- 
native holds. 
Either 



i) the system (1.2) is equivalent to the system (1.51 



ii) there exists a polynomial solution to 

div(AVw 1 + BVu 2 ) =0 , m B, 
div(CVw 1 + DVu 2 ) = , on B, 

which is a homeomorphism of a closed disk B onto u(B) and such that in the center 
of the disk, denoted by O, we have 

det Du(O) = 0. 

Remark 1.7. It may be evident that the Rado-Kneser-Choquet and the H. Lewy 



Theorems continue to hold for any system of the form ( 1.5 1, since it can be elemen 



tarily reduced to a Laplacian diagonal system via a linear change of the independent 



coordinates. It is although rather remarkable that Theorems 1.4 1.6 show that the 
Rado-Kneser-Choquet and the H. Lewy Theorems do not extend to very slight 
perturbations of the Laplacian diagonal system such as, for instance, the following 
one 

u xx u yy 0, 

(1 + e)«L + u lv = °> 
where e is any positive number. 
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2. Proofs 



In what follows, when no ambiguity occurs, we shall identify points (x,y) £ 
x 



with column vectors 



y 



Also, for 9 € K, we shall denote eg = cos 9 , sg = sin 9. 



For the proofs of Theorems |1.4| and |1.6| we shall make use of the following two 
propositions, which we will prove at the end of this Section. 



Proposition 2.1. Let the ellipticity condition (1.3) be satisfied. If the system (1.2) 
is not equivalent to (1.5), then there exists 9 € [0,27r] and a quadratic polynomial 
p(x, y) — \(ax 2 + 2bxy + cy 2 ) such that 



co 
■so 



se 
co 



x* + y* 

p(x,y) 



is a solution to (1.2). 



Proposition 2.2. Let the ellipticity condition (1.3) be satisfied. If the system (1.2) 
is not equivalent to (1.5), then there exists 9 € [0,27r] and a cubic polynomial 



q{x,y) 



+ bx 2 y + cxy 2 + d 



y 



such that 



( u 1 N 


)=( 


' eg 


~: ) ( 


1 u 2 , 




v so 





x(x 2 + y 2 ) 



is a solution to (1.2). 



Proof of Theorem \1.4\ We assume that (1.2 1 is not equivalent to (1.5 1. We choose 
9 according to Proposition 2.1 and we write 



co 
*o 



~Sg 
CO 



Being linear mappings solutions to ( 1.2 ), we have that also the following is solution 
to ( |0| 

x 2 + y 2 — 1 
ky + p(x,y) 

where k ^ is a constant to be determined later on. 



Ro 



We choose B 



{(x,y) £ M 2 : (x-\f + y 2 < |}. We have 
x = -l + x 2 + y 2 on dB. 



We now select <&. We set 



M = 



1 
k 



V(x,y) = 



y+ lp(x,y) 



&(x,y) = R e M^(x,y). 



A straightforward calculation shows that when 
(2.1) fc>|6| 

^ is a homcomorphism of dB onto a closed convex curve T. Consequently <i> is also 
a homcomorphism of dB onto the closed convex curve 7 = RgMT. 
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Let D be the bounded convex domain such that dD = 7. It is easy to check that 
DC R e MS where 



However 



S = i (v\v z ) € R : — < < 



»(().()) = 1/ ( J ) 7? /?„.US'. 



□ 



Remark 2.3. Note that we can compute the the Jacobian determinant of u in terms 
of the coefficients of p and obtain det Du = 2kx + 2(b(x 2 — y 2 ) + (c — a)xy). Hence 
the Jacobian determinant vanishes at (0, 0) and, in fact, it changes sign across its 
nodal line 

H = {(x, y) G M 2 : kx + (b(x 2 - y 2 ) + (c - a)xy) = 0} . 

which is always an hyperbola (unless a = b = c = when it degenerates in a straight 
line). Note that kx + (b(x 2 — y 2 ) + (c — a)xy) is positive in (1, 0) zero in (0, 0) and 



negative in (—1,0) because of (2.1|. See Figures 1 and 2 where the behaviour of u 



and its Jacobian determinant are depicted in the specific case of the Lame system. 



Proof of Theorem \1.6\ We assume again that (1.2 1 is not equivalent to (1.5 1. We 



choose 9 and q according to Proposition 2.2 Being linear mappings solutions to 



/ u 1 x 


)=( 




7 ) ( 


{ u 2 , 









(1.2 1, we have that also the following is solution to (1.2 1 

x(x 2 + y 2 ) 
y + l{x,y) 

It is now easy to check that, once we have chosen the polynomial q according to 
Proposition |2.2| the following two properties hold. First there exists a positive 
radius r such that one has det Du(x,y) > if (x,y) G B r (0)\{0}, where we have 
set O = (0, 0), and second det Du(O) — 0. Choose < p < r, and denote by 

* = u \dB p (oy 

A very simple calculation shows that <!> maps dB in a one to one way onto a closed 
curve 7 provided one has 

2 + bp 2 > and 3 + dp 2 > 0. 

Let D be the bounded domain such that dD = 7. We can now apply a topological 
result of Meisters and Olcch [20 . Indeed we have a smooth mapping u defined on 
a closed disk B and which is a local homeomorphism at each point of B, with the 
possible exception of the point O only. Moreover, the restriction of u to dB is a 
homeomorphism. Therefore the hypotheses of Theorem 1 in [20 (see also Corollary 
2) are satisfied and we can conclude that u is a homeomorphism of B p {0) onto 
D. □ 

Proof of Theorem\l.S\ It suffices to verify that the Lame system is not equivalent 



to any elliptic system of the form ( 1.5 ). In fact it can be rewritten in the form ( 1.2 ) 
with the following choices 
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Now B, C and D can be scalar multiples of A only if fj, + A = 0, which contradicts 
the assumption fj, + A > 0. □ 

Remark 2.4. Note that our assumptions /i>0,/i + A>0 correspond the the 
strong convexity assumption, which, as is well known, is stronger that the ellipticity 



condition ( 1.3 ) 



In the following picture we illustrate our results for the case of the Lame system. 



We choose A = /z = 1. With this choice, condition (2.1l takes the form k > 
(l+-\ZlO) and, for the picture we have chosen the limiting value k = 2(l + vT0). 
Note that, in order to facilitate visibility, the coordinates in the u 1 and the u 2 
directions are scaled differently. 




Figure 1. dB and its image <&(dB). 




Figure 2. Left: circles C r of varying radii and the nodal line of 
the Jacobian (an hyperbola) drawn within B. Right: the images 
U(C r ). 



Proof of Proposition \2.1\ We look for 6 S [0, 27r] and a quadratic polynomial p(x, y) 
\{ax 2 + 2bxy + cy 2 ) such that 

,i 



(2.2) 



ce -se 

S0 Cg 



x 2 + y 2 



is a solution to (1.2). 

Since the Hessian matrices of u l are constant, the system (1.2 1 is equivalent to 
the following two equations 



( -s,,F - c„G) | b 



Y 
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where 



and Y = 



F = 

Y 1 
Y 2 



an 
en 



2a i2 
2ci 2 



ai2 
C22 



G 



hi 
da 



2b 12 
2d 12 



b'22 

d 22 



is a vector of known data, possibly depending on 9. Given 9, 



the above system has at least one solution if the rank of —sgF + cgG is two and 



consequently we find a solution of the form (2.2 1 to ( 1.2 1. If this were not the case, 
then for every 8 there exists 

(2.3) 

Choosing 9 



such that 

C<t,(-SeA + CgB) + S<j,(-SgC + CgD) = 0. 

yields that B and D are linearly dependent. Similarly, by choosing 



6 = w, we get that A and C are linearly dependent. Recalling that, by ellipticity 
(1.3 1 j A and D are positive definite, and hence nontrivial, we obtain that B = o~D 
and C = j A for suitable constants a, 7 € R. Plugging these linear dependencies 
into (2.3 1 and using once more ellipticity, we obtain that also D is a scalar multiple 



of A. In conclusion B, C and D are scalar multiples of A. 



a 



Proof of Proposition 2.2 We look for 9 € [0, 27r] and a cubic polynomial 

„3 



such that 



/ u 1 N 


H 




7 ) ( 


I « 2 , 




v s 





+ bx 2 y + cxy 2 + d^— 



x(x 2 + y 2 ) 



<l{x,y) 



is a solution to ( 1.2 1. This time the Hessian matrices of u l are of the form xHi+yH 2 



for suitable constant matrices Hi and H 2 and thus ( 1.2 1 imposes the following four 
conditions. 

c e G) 



{s e F- 



b 



= Y 



where 



F = 





an 


2oi2 


022 


\ 







an 


2a i2 


«22 




Cll 


2C12 


C22 





V 





Cll 


2C12 


C22 / 



G = 



1 bu 


2&12 


b 22 


° \ 





fen 


2bi 2 


^22 


dn 


2di2 


d 22 





V 


dn 


2di 2 


C^22 / 



and Y £ R 4 is a data vector. We make use of the following linear algebra fact. 
Given any two 2x2 symmetric matrices 



M = 



run 
m i2 



mi 2 
m 22 



and S = 



Sll 

S12 



if we have 



det 



/ 


mn 


2mi2 


m 22 


\ 







mn 


2mi2 


m 22 




Sll 


2S12 


s 22 





V 





Sll 


2S12 


S22 / 



S12 
S22 



0, 



and either M or S is positive definite, then M and S are linearly dependent. This 
fact may be verified in many ways, for instance with the aid of Gaussian elimination. 



MATERIAL INTERPENETRATION 



9 



Assume that (—sgF+cgG) is singular for every 9, then we must have that the matrix 
(— sgC + cgD) is a scalar multiple of (— sgA + cgB) at least for all those 9 for which 
{—sg A + cgB) is positive definite. Equivalently, switching the roles of A,B and 
D, C, if (—sgC + cgD) is positive definite, then (— sgA + cgB) is a scalar multiple of 
(—sgC + cgD). Recalling that the matrices A, D are positive definite by ellipticity, 
and choosing 9 = 0, we deduce B — oD and C — jA for suitable constants 
cr,7 € i. Moreover, it is evident that there exists an open interval / containing 
^ for which (— sgA + cgB) remains positive definite as long as 9 E I. Thus, for 
all 9 E I, there exists 4> such that (2.3 1 holds true, and from now we can argue 



similarly as in the proof of Proposition 2.1 □ 
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